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I. INTRODUCTION 

This paper is a continuation of the work whose 
focus was on the simplest characteristics (level density) 
of a weighted sum of unitary random matrices — a non- 
Hermitian random matrix model encountered e.g. in 
quantum information theory or the theory of random 
walks on trees. Inspired by the first of these areas 
of physics (Sec. II CI) . I introduce more general models 
(Sec. II Ap and calculate — using the tools of free probabil- 
ity (Sec. II B| ) — their level densities (Sec. HH) . 



A. Models 

1. Definitions of the models 

The goal of this article is a basic study of the following 
non-Hermitian random matrix models T, W, V, which 
I will call "generalized Bures products" for the reason 
explained at the end of Sec. II C 21 In order to introduce 
them, I recall two other matrix models, S and P, inves- 
tigated recently in the literature: 

(i) Let L > 2 be an integer; let U;, for I — 1,2, . . . , L, 
be independent N x N unitary random matrices belong- 
ing to the Wigner-Dyson "circular unitary ensemble" 
(CUE or "Haar measure"; i.e., their eigenvalues are on 
average uniformly distributed on the unit circle); let Wi 
be arbitrary complex parameters. Then define a weighted 
sum, 



S = wiUi + W2U2 



(1) 



This non-Hermitian model may be referred to as the 
"generalized Kesten model" (cf. also 0,11), and has 
been considered in some detail in [ij. 

(ii) Let K > 1 be an integer; let A^, for 
k — 1,2,...,K, be rectangular (of some dimensions 
Nk X Nk+i) complex random matrices such that all the 
real and imaginary parts of the entries of each ma- 
trix are independent random numbers with the Gaus- 
sian distribution of zero mean and variance denoted by 
al/{2{NkNk+i)^^^), i.e., the joint probability density 
function (JPDF), 



K 

P{{Ak}) cx []exp 

fc=i 



Tr(AtAfc)); (2) 



these are rectangular "Ginibre unitary ensembles" 
(GinUE) 0, Q. Then define their product, an 
Ni X Nx+i random matrix, 

P = AiA2...A^f. (3) 

This model has been thoroughly investigated (cf. e.g. 0- 

Ei). 

Now, having a number of copies of S and P (all of them 
assumed statistically independent from each other), one 
may use them as building blocks of various products: 

(i) A product of J > 1 weighted sums ^ , 



T = S1S2...S,, 



(4) 



where each sum Sj has arbitrary length Lj and complex 
weights Wji. 

(u) A product of (g]) and 



W = TP, 



(5) 
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i.e., this is a string of J generalized Kesten ensembles 
and K Ginibre unitary ensembles. Note that one has to 
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set N — Ni, and then W is rectangular of dimensions 
iVi X Nk+1. 

(iii) Finally, a most general string of generalized Kesten 
and Ginibre unitary ensembles is obtained by multiplying 
a number / > 1 of random matrices ([S]), 



V = W1W2 



. w 



(6) 



The dimensions of the terms {Ni^i x Ni^Ki+i, for i — 
1,2, ... ,1) must obey Ni^Ki+i = ^i+i.i, and then V is 
rectangular of dimensions A^i^i x Nj^Ki+i- 



2. Thermodynamic limit 

The techniques I apply to approach the above models 
are valid only in the "thermodynamic limit," i.e., for all 
the matrix dimensions infinite and their "rectangularity 
ratios" finite, 



N^N,,N2, 



Nk- 



. .,Nk+i -> 00, 
finite. 



(7) 



[I use here the notation for the model W (O; for V 
one should add an index i.] 



3. Mean densities of the eigenvalues and singular values 

Level densities. This paper is just a most basic study of 
the models X e {T, W, V}, namely, I will be calculating 
only: 

(i) The mean density of the eigenvalues ("mean spec- 
tral density" or "level density"). 



1 ^ 

p^{z,z)^-Y,{5^^^ (z-A.)), 
1=1 



(8) 



where the mean values are with respect to the JPDF of 
X and are evaluated at the complex Dirac delta functions 
at the (generically complex) eigenvalues \i of X. Note 
that one must set iV^'+i = N (i.e., ri 1) for W to be 
a square matrix. 

(ii) The mean density of the "singular values," which 
are the (real and non-negative) eigenvalues /i^ of the 
Nk+1 X Nk+i Hermitian random matrix H = X^X, 



P-Si{x) 



1 



Nk + 1 



(9) 



where the Dirac delta is real. This time, ri may acquire 
any positive value. 

Future work — universal quantities. Certainly, this is 
but the first step in understanding the considered mod- 
els. The level density of a random matrix is known not 
to be universal, i.e., it depends on the precise probabil- 
ity distribution of the matrix. (However, cf. the end of 



Sec. II B 3l for a hint in favor of a certain universality of the 
level densities of our models.) As a next step, it would 
be desirable to investigate some "universal" properties 
of our models, i.e., depending solely on their symmetries 
but not the specific probability distributions. Some basic 
universal observables would be a "two-point connected 
correlation function," 



connected 



{x,y) = piiix,y) - pii{x)pii{y), 



(10a) 



- ■ n -|- J- 

piiix,y) = ^ — E - fJ-r) S {y ~ Pj)) 

(10b) 
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(here written in the Hermitian case), or the behavior of 
the level density/correlation function close to the border- 
line (for non-Hermitian models) or the edges (for Hermi- 
tian models) of the spectrum. 

Universal erfc scaling. Even though I focus on the level 
densities in the bulk of the spectrum and not close to its 
borders, the models T, W, V share a certain property 
(their mean spectrum is rotationally-symmetric around 
zero; cf. Sec. II B 2[) which permits an application of the 
so-called "erfc conjecture" (cf. Sec. lIB 6)) . which describes 
the universal way in which the mean spectral density 
is modified close to the borderline of the "mean spectral 
domain" V (i.e., the subset of the complex plane where 
the eigenvalues of an infinite random matrix fall). The 
pertinent form-factor psp depends on one or two param- 
eters, and I positively test this hypothesis on a number 
of examples for each class of the models T, W, V by 
fitting these parameters to match the Monte Carlo data. 
However, it still remains to prove this hypothesis, as well 
as analytically determine the form-factor parameters. 

Future work — microscopic quantities. Another step in 
the research on our models could be to consider finite 
matrix dimensions, and to calculate the complete JPDF 
of the eigenvalues and singular values of S, P, T, W, V. 
This however is a much more involved task. 



B. Tools 

In this Section, I sketch the means by which the level 
densities ^] of the models [(U), (O, dH)] wih be 

evaluated in the thermodynamic limit ([T])- I do not delve 
into details, as they have been addressed e.g. in (ll. [l2j. 



1. Basic language of random matrix theory 

Mean spectral densities and Green functions. First of 
all, it is convenient to replace the mean spectral density 
of any {N x N, N 00) Hermitian random matrix H 
or non-Hermitian random matrix X with an equivalent 
but more tractable object in the following way: Since the 
definitions [Q, (HI)] exploit the real (Hermitian case) or 
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complex (non-Hcrmitian case) Dirac delta function, one 
uses their respective representations, 



5{x) = —- — r lini ( : :- 

zvri e^a+ \ x + i€ X — le 



1 



(11a) 

(lib) 



which prompt one to introduce the following "holo- 
morphic" or "nonholomorphic Green functions (resol- 
vents)" [lM3, 



1 



1 ^ 



.iTr 

N 



1 



zl 



N 



(12a) 



1 ^ 

Gx(z,z) = lim lim — > 

1=1 

= lim lim — Tr- 

e-i.OAr-i.oo N 



z - Xi 



\z-Xr + e^ 



zIn - Xt 



(zlAr-X)(zlAr-Xt)+e2lAr 



(12b) 



Thanks to (|lla|) - (jllbll . the densities are straightforward 
to reproduce from them. 



Pii(x) 



— lim (GH(a: + ie) - GH(a; - ie)) . 

ZTTl e-i.O+ 



P:x.iz,z) = -chGyi{z,z) 

TT 



(13a) 
(13b) 



while the resolvents prove to be handier from the com- 
putational point of view. 

Remark: Eqs. [ (fT2b| . (pb]) ] are valid for z inside the 
mean spectral domain T). Outside it, the regulator e 
may be safely set to zero, and the nonholomorphic Green 
function reduces to its holomorphic counterpart (I12al) . 



Gx(2,z) = Gx(^), ioiz^V. 



(14) 



An implication is that by calculating the Green function 
both inside 2? and outside T), and then equating them 
according to (|14p . one arrives at an equation of the bor- 
derline dT>. 

M -transforms. In addition to the resolvents, one of- 
ten finds even more convenient to work with their simple 
modification, the "holomorphic" or "nonholomorphic M- 
transforms " 



Mh(z) = zGu{z) - 1, 
M^(z,z) EE zGyi{z,z)-\. 



(15a) 
(15b) 



The Hermitian M-transform (|15ap has the interpretation 
of the generating function of the "moments" toh.m (if 
they exist), by expanding around z = oo, 



Mu{z) 



TV 



Tr (H") . (16) 



2. Rotational symmetry and N -transform conjecture 

Rotational symmetry of the mean spectrum. In this 
paper, I investigate only non-Hermitian random matrix 
models X with a feature that their mean spectrum is 
rotationally-symmetric around zero, i.e., that their den- 
sity ([U depends only on 



R^\z\. 



(17) 



Equivalently, this means that their nonholomorphic M- 
transform (Il5b[) also depends only on the radius. 



Mx(z,z) =9[Hx (i?' 



(18) 



giving the relevant radial part of the density through 



p5t^ (i?) = 2^i? p^{z,z) 



(19) 



N -transform conjecture. For such ensembles, it has 
been suggested and numerically confirmed on a number 
of examples [H, [HI, [l3] that there exists a simple rela- 
tionship between the mean densities of their eigenvalues 
and singular values (i.e., eigenvalues of X^X) — the "A^- 
transform conjecture." In order to express it, one needs 
to define the functional inverse of the M-transform. In 
the Hermitian case (|15ap . it can be directly done. 



AfH(^H(^)) = Ah(Mh(^)) 



(20) 



obtaining the "holomorphic A^-transform." In a generic 
non-Hermitian case (jl5bp . however, such an inversion is 
obviously impossible. But with the rotational symme- 
try present (1181) . the M-transform again depends on one 
argument, and the inversion becomes generically doable. 



5Hx (OTx(^)) 



OTx (OTx (i?')) = R\ (21) 



which is called the "rotationally-symmetric nonholomor- 
phic A^-transform." [To be more precise, ?tx(-z) defined 
by the left equation in (j21l) is a holomorphic continu- 
ation of the one defined by the right equation.] Now 
the hypothesis states that these A^-transforms of X and 
H = X^X remain in a simple relation. 



^xtxW 



z + l 



^x(^). 



(22) 



Typically, one of these random matrices will be much 
more accessible by analytical methods than the other, 
and therefore, Eq. ([^^ will provide a way to that more 
complicated model. 

I will henceforth make extensive use of this hypothesis; 
numerical tests of the results thus obtained will further 
support it. 



Free pro 



and model P 



Algorithm for computing the densities o/P and 'P'^'P . 
The model P ([3]) is a product of independent (rectangu- 
lar) random matrices A^. If ri ^ 1 (i.e., P is rectangu- 
lar), then one is interested in the eigenvalues of P^P only. 
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while if ri = 1 (i.e., P is square), then its mean spectrum 
has rotational symmetry around zero (O 53, El, and 
can therefore be related (|22p to the spectrum of P^P. 
The level density of this latter matrix can now be ex- 
pressed (by cyclic shifts of the terms) [TTJ through the 
density of the product of the Hermitian ("Wishart" |18|]) 
matrices Aj^A^. Hence, the situation is suitable for the 
employment of the "multiplication law" of free probabil- 
ity calculus of Voiculescu and Speicher [19, 20,]. 

Multiplication law in free probability. This theory is es- 
sentially a probability theory of noncommuting objects, 
such as large random matrices. Its foundational notion is 
that of "freeness," being a proper generalization of sta- 
tistical independence. Qualitatively speaking, freeness 
requires the matrices Hi, H2 to not only be indepen- 
dent statistically, but also "rotationally" (i.e., no distin- 
guished direction in their probability measures, i.e., de- 
pendence only on the eigenvalues) — independent random 
rotations from the CUE, UiHiUj and U2H2U2, would 
ensure that condition if necessary. 

Assuming freeness, there is a very useful prescription 
for computing the mean spectral density of the product 
H1H2 (assuming it is still Hermitian, which obviously 
is not always the case; but the same procedure applies 
when all the matrices are unitary), once the densities of 
the terms are known: (i) Encode the densities by the 
holomorphic Af -transforms (|15ap . Mhi(z) and Mh2(^), 
as outlined in Sec. II B II (ii) Invert them functionally 
to find the respective A^-transforms ((20|) . iVHi(z) and 
^H2(-^)- (iii) The A^-transform of the product follows 
from the "multiplication law," 



HiH 



z + 1 



(23) 



(iv) Invert the result functionally to find the M- 
transform, and thus the mean spectral density of the 
product. 

Eigenvalues and singular values o/P. Implementing 
the above steps to P, one readily discovers 



7Vptp(z) 



a^^AT^ n ( - 



where for short, a = HfcLi '^fc, from which stems both the 
mean density of the singular values [^U^, (|15ap . (|13al) ] 
and (provided that ri = 1) the eigenvalues [ d^^ . (I?!]) . 
P^ ] of P [in the latter case, the mean spectral domain 
I? is a centered disk of radius Rp.^t. = a, cf. P7al) - (j37bp ]. 

Furthermore, the authors of conjecture and numer- 
ically affirm that (j24p is to some degree universal, namely 
that it holds for the matrix elements of the A^ not only 
IID Gaussian but also arbitrary IID obeying the Pastur- 
Lindeberg condition. Depending on a possible analogous 
universality of the level density of S, also the densities of 
T, W, V may exhibit a certain universality. 

Let me close by saying that an algorithm parallel to 
the one described above will also be applied to derive the 
main results of this paper (cf. Sec. Ill A| ). 



z 



(24) 



4- Quaternion free probability and model S 

Although the topic of this work is multiplication rather 
than summation of random matrices, I should mention 
another [in addition to the technique outlined in Sec. II B 3l 
and Eq. (IM)) ] pillar of this article — the results for the 
model S ^ derived in 

Addition law in free probability. In the realm of Hermi- 
tian random matrices, free probability provides — besides 
the multiplication law — a rule for summing free ma- 
trices Hi, H2. This time, one should invert functionally 
not the holomorphic M-transforms, but the Green func- 
tions 

Gh(Sh(z)) =Sh(Gh(2)) -z, (25) 

which is known as the "holomorphic Blue function" 12 il l 
This objects then satisfies the "addition law" [l^ [Ic 



Shi+H2(2) = Shi(z) 



z 



(26) 



upon which it remains to functionally invert the left-hand 
side, which leads to the holomorphic Green function of 
the sum, and consequently, its mean spectral density. 
(Recall that in classical probability theory, an analogous 
algorithm makes use of the logarithms of the character- 
istic functions of independent random variables.) 

Addition law in quaternion free probability. For non- 
Hermitian random matrices, a construction parallel to 
(I25|)-(l26l) has been worked out in [H H^. Basically, it 
consists of three steps: 

(i) "Duplication" [H (cf. also [H, [13)— the non- 
holomorphic Green function (I12bp has a denominator 
quadratic in X, which makes its evaluation hard, and one 
needs to linearize it by introducing the (2 x 2) "matrix- 
valued Green function," 



Qx(z.z) = lim lim — bTr/— —3 — r 

^ ' ^ e^O N^oo N \ ® 1 AT - Xdupl- 



where for short. 



z le 
ie z 



■^dupl. 



X On 
On Xt 



and the "block-trace," 



bTr 





B ^ 


l^f TrA 


TrB \ 




D J 


' - I TrC 


TrD ) 



(27) 



(28) 



(29) 



In other words, this object resembles the holomorphic 
Green function p2ap . and thus can be approached 
by methods designed for Hermitian random matrices; 
the cost being the need to work with 2x2 matri- 
ces rather than complex numbers. The nonholomor- 
phic Green function lies precisely on its upper left entry, 
[^/x(z,z)]ii = Gx(z,z). 

(ii) The extension to the complete quaternion space — 
by replacing in (j27p the infinitesimal e with a finite com- 
plex number. 



c id 
id c 



c, d G C, 



(30) 
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one obtains a quaternion argument of the 
Green function," 



1 



-bTr 



1 



Xdupl. 



'quaternion 



(31) 



The mean spectral density is reproduced from this 
quaternion function by approaching the complex plane 
(c — z, d — e), just as in the Herniitian case, it is ob- 
tained from the complex Green function by approaching 
the real line {z = x ± ie) p3ap . 

(iii) The functional inversion can now be performed in 
the quaternion space, 



ax(Sx(Q)) = Sx(ex(Q)) = Q 



(32) 



(the "quaternion Blue function"), and the "quaternion 
addition law" for free non-Hermitian random matrices 
can be proven. 



(33) 



6xi+x.(Q) =Sxi(Q) + 6x.(Q)- Q"^ 



Eigenvalues and singular values of S. Formula (j33p 
massively simplifies calculations of the mean spectral 
density of sums of free non-Hermitian ensembles. In par- 
ticular, harnessing it to the weighted sum ([T]) implies 
that: (i) Its mean spectrum exhibits the rotational sym- 
metry around zero (jl8l) . (ii) Its rotationally-symmetric 
nonholomorphic A^-transform (PT|) is a solution to the 
following set of {L + 2) polynomial equations. 



models with the rotational symmetry (|18p present, their 
mean spectral domain V is always a disk or an annulus. 
It has been proven for the models S and P; for T, W, 
V, I will assume it holds, and support it a posteriori by 
numerical simulations. 

Radii of the disk/annulus. If this conjecture is true, 
then the external and internal radii of the annulus (be- 
coming a disk if the internal radius shrinks to zero) are 
still given by (I35ap - (|35bp . provided that there are no 
zero modes. Indeed, on the borderline of the domain 
T), the nonholomorphic Green function (from the inside 
of V) and the holomorphic one (from the outside of V) 
must match But the holomorphic Green function 

compatible with the rotational symmetry (flSl) has the 
form Gx(2;) = (1 + 97tx(-R^))/2, hence, holomorphicity 
leaves one possibility, i.e., that S[Hx(-R^) is a constant, 
G'x(z) — (1 + ^)/z. In the external outside of T) (in- 
cluding z ~ cx)), it is well-known that G'x(z) ^ 1/z as 
z — > oo, which means that 93t = 0. In the internal outside 
(including z = 0), supposing there are no zero modes, 
the Green function must be just zero, i.e., 3Jl = — 1. If 
there are zero modes, Px™ ™°'^'''^(z,z) = a5'^'^\z ^'z) , they 
are obtained by applying (jl3bp to the Green function 
Gf ° m°d<=«(z) ^a/z regularized according to Eq. (fTTb| : 
therefore, 93T = a — 1. In other words, 



R 



ext. 
2 



^x(O), 
OTx(a- 



!)• 



(37a) 
(37b) 



1=1 



-C = 



z{z + l) 
OTs(^) ' 
Ml {Ml + 1) 



/ = 1,2,...,L, 



(34a) 

(34b) 
(34c) 



where G > and Mi are auxiliary unknowns, (iii) Its 
mean spectral domain 2? is either a disk or an annulus, 
whose external and internal radii are found from 



Rt 



9ts(0), 



=9Ts(-l). 



(35a) 
(35b) 



Furthermore, Eq. (|22p means that the singular values are 
derived from an identical set of equations except (j34bp 
which turns into 



-G 



^sts(^)' 



(36) 



e ring conjecture 



Single ring conjecture. It is another hypothesis [l| (be- 
ing a generalization of the "Feinberg-Zee single ring the- 
orem" (27l - [30| ) that for non-Hermitian random matrix 



6. erfc conjecture 

Calculations in the thermodynamic limit ([7]) are ca- 
pable of reproducing the mean spectral density only in 
the bulk of the domain T>, but not close to its border- 
line. However, based on earlier works [1, |3TI - [33| . it has 
been suggested [l|, [111, [131 ^^"^ extend the radial mean 
spectral density (fT9|) to the vicinity of the borderline, 
provided that the rotational symmetry (|18p holds. 



For each circle R — Rh enclosing T) (one or two; 
cf. Sec. II B 5)) . one should simply multiply the radial den- 
sity Px"^ (i?) by the universal form-factor. 



/jv,gb,_Rb,sb(i?) = ^erfc (qhSh {R ~ Rh) 



(38) 



where erfc(a;) = J^dtexp(— t^) is the complemen- 
tary error function, while the sign st, is -|-1 for the ex- 
ternal borderline and —1 for the internal borderline, and 
(7b is a parameter dependent on the particular model, 
whose evaluation requires truly finite-size methods, but 
whose value I will adjust by fitting to the Monte Carlo 
data. 
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C. Motivation 

1. Interesting mathematical properties 

I find the models T, W, V mathematically inter- 
esting for the following reasons: (i) They are non- 
Hermitian, and the theory of such random matrices is 
richer and much less developed than for the Hermitian 
ones (cf. e.g. [3l|). (ii) They belong to a special class 
of non-Hermitian matrices, namely, with rotationally- 
symmetric mean spectrum (jl8p . As such, they con- 
jecturally exhibit certain features, which demand test- 
ing (and eventually proofs): the iV-transform conjec- 
ture (Sec. II B 2p . which makes them reducible (at least 
concerning the level density) to Hermitian models; the 
single ring conjecture (Sec. IIB 5|) : the erfc conjecture 
fSec. IIB~6l) . (iii) They contain two operations widely in- 
vestigated in the literature on random matrices, albeit 
rather separately: summation (cf. e.g . | M [1 , [22| - [27l Is^- 
[lil) and multiplication (cf. e.g. [Mlllilll). 



2. Applications to quantum information theory 

The models T and W are encountered in the theory 
of quantum information. This application has been de- 
scribed in m, HO], as well as in Sec. I B 2 of and 
I refer the reader to these sources for a more detailed 
exposition (the textbook ,51j is an introduction to the 
subject). It is debatable whether one can find the model 
V in this theory; if not, the reader may treat it just as a 
mathematically natural generalization of W. 

Density matrix. Fundamental objects in quantum in- 
formation theory are "mixed states," i.e., statistical en- 
sembles of quantum states. Such random states arise in 
a number of important settings: (i) If a quantum system 
interacts with its environment in a complicated (noisy) 
way, which may be regarded as random, (ii) If a quantum 
system is in thermal equilibrium, (iii) If the preparation 
history of a quantum system is unknown or uncertain 
(such as for quantum analogues of classically chaotic sys- 
tems), (iv) If one investigates generic properties of an 
unknown complicated quantum system, one may assume 
it is random, (v) If a system consists of subsystems which 
are entangled, each of them must be described by a mixed 
state (and quantum entanglement is a central feature in 
the theory of quantum computers). A classic example 
is light polarization: a polarized photon can be written 
as a superposition of two helicities, right and left circu- 
lar polarizations, {a\R) + b\L)) (a "pure state"); whereas 
unpolarized light may be described as being \R) or |L), 
each with probability 1/2 (a mixed state). 

A mixed state cannot be represented by a single state 
vector — a proper formalism is that of a "density matrix." 
If one considers a statistical mixture of N pure states 
each with probability pi G [0, 1| {J^iLiPi = 1): the 
density matrix is defined a.s p = J2i=iPi\''Pi){'^i\ 
vex combination of pure states; called sometimes an "in- 



coherent superposition"), and it has this genuine prop- 
erty that the expectation value of any observable A is 
given by Tr(/3A). More generally, any operator which is 
Hermitian, positive-semidefinite (its eigenvalues are non- 
negative) and has trace one (its eigenvalues sum to one) 
may be considered a density matrix. 

I will be interested in complicated composite quantum 
systems. For example, for a bi-partite system consisting 
of a subsystem A of size A^i and B of size N2 (with or- 
thonormal bases, {|j)^} and {|j)b}), a general pure state 
of the full system is "entangled," i.e., it cannot be written 
as a tensor product of pure states of A and B, 



Ni N2 



(39) 



1=1 j=i 



In other words, Aov B cannot be said to be in any definite 
pure state. But they can be characterized by a density 
matrix — consider any observable A on A; its expectation 
value in the state reads Tr(p_^A), where the "reduced 
density matrix" of A, 



Pa 



xxt 



Tr(XXt)' 



(40) 



and X is rectangular iVi x N2 complex matrix. (This 
may equivalently be obtained by performing the "par- 
tial trace" Trg of the density matrix of the full sys- 
tem, i-e., J2j B{jH){ip\j)B, and normalizing it.) 
Therefore, the eigenvalues of (which are the properly 
normalized singular values of X) are precisely the above 
probabilities pi. This picture should be supplied by the 
discussed above basic notion of replacing the complicat- 
edness of the full system by randomness — considering X 
to be some random matrix, and calculating the mean 
density of its singular values. 

An important measure for a mixed state is its "von 
Neumann entropy," 



JV 



S{p) = -Tr(plogp) = - ^pilogpi 



(41) 



It represents the degree of randomness (mixture) in the 
mixed state (thus, quantum information; thus, it also 
measures entanglement) — it is the larger, the more dis- 
perse the probabilities pi are; it is zero for a pure state, 
and reaches its maximum logA^ for all the probabilities 
equal [the full system is in the "maximally entangled 
(Bell) state"]. Moreover, a measurement can never de- 
crease this entropy; consequently, a measurement may 
take a pure state to a mixed one, but not conversely 
(except there is a greater growth of entropy in the envi- 
ronment). 

I will focus on some physically motivated "structured 
ensembles" of random states in which an important role 
is played by the tensor product structure of the Hilbert 
spaces of the subsystems, i.e., which are invariant under 
local unitary transformations in the respective Hilbert 
spaces. 
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Model S. The following construction leads to the ap- 
pearance of the model S in the above setting: 

(i) Consider a bi-partite system {A, B) of size N x N 
in the Bell state, 



1 ^ 



(42) 



(ii) Consider the following type of randomness — apply 
to this Bell state L independent random local unitary 
transformations (belonging to the CUE) U; in the prin- 
cipal system A, and form a coherent superposition of the 
resulting (maximally entangled) states with some weights 
wi, 



1=1 



(S®ljv)|vl/+^) 



(43) 



1 ^ 



10-4 «) 



(iii) Hence, the reduced density matrix for A is given 
by (|40)) with X = S. [The normalization constant is 
Tr(SSl') = l^'P + • • where the dots are much 

smaller than N in the large- iV limit ([7]). So one may 
set for convenience, X]i=i l"*^'!^ — li ^'^d investigate the 
singular values of S rescaled by A^.] 

Model T. It is a matter of a nested repetition of the 
above procedure to obtain the model T. For instance, 
for J = 2: 

(i) Consider the above pure state , constructed ac- 
cording to a matrix S2 of length L2, and take its Li 
copies. 

(ii) Perform a random independent CUE rotation Ui; 
of each of these copies, and form their coherent superpo- 
sition with some weights wu . 

(iii) The reduced density matrix for A PO]) will have 
X = T = S1S2. One may continue this process any J 
times. 

Model P. This ensemble originates from a measure- 
ment in the product basis of Bell states, in the following 
way: 

(i) Consider a composite system with 2K subsystems 
of sizes 

Al,^ A2, A3^, ■ . . , A2K-2,A2K-U A^ (44) 
size A^i size A^2 size Nk ^i^^ Nk+\ 

(ii) Take an arbitrary product state 
1^0) = |0)^, (55 |0)^, . . . ® |0) 

A2K J and apply to it 
random unitary local transformations acting on the 
pairs of subsystems. 



2K-1 K jIV'o). (45) 



The result is separable with respect to the above pairing, 
i.e., it can be expanded in the product basis as 



Ni N2 

E E • 



Nk Nk+1 

E E 



■ VllAi ® \i2)A2 ® l«2)-43 (5 



(46) 



i...(E) \iK+i)A2K, 



where the coefficients, gathered into K rectangular 
{Nk X Nk+i) matrices A^, may be assumed independent 
Gaussian 

(iii) Consider the Bell states p2)) on the pairs (^^2, ^3), 
• • . , iA2K-2,A2K-i), and project |-0) onto their product. 



K 



r = i. 



\k=2 



(47) 

which leads to a random pure state describing the re- 
maining two subsystems, Ai and A2K, through the ma- 
trix P, 



oc 



Ni Nk+1 

E E 



JjriK+ihl)-4i ® Vk+i)A2k- 



(48) 



ii — 1 ii. 



(iv) The reduced density matrix for Ai (i.e., the nor- 
malized partial trace over A2K) is therefore with 
X = P. 

Model W. A direct combination of the above two al- 
gorithms leads to the model W — one should consider the 
random pure states corresponding to the model T on the 
2K subsystems ((44)) . and proceed as above. 

One may find in the literature only an expression for 
the mean density of the singular values of W in the case 
of J = 1, Li = 2, wii = I/V2 (for / = )-^)^ = 1' 
CTi = 1, ri = 1 — the "Bures distribution" l52l. l53|. 




(49) 



for X € [0, /3] and zero otherwise, where for short, 
/3 = 3^3. The chief purpose of this paper is to extend 
this result to arbitrary values of the parameters, even to 
a product of a number of matrices W, as well as to the 
mean density of the eigenvalues. 



II. GENERALIZED BURES PRODUCTS 

A. Generalized multiplication laws 

This paper deals with products of random rectangu- 
lar or non-Hermitian matrices, hence, I will start from 
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adjusting the free probability multiplication law (I23p to 
such a situation. Consider a product of arbitrary inde- 
pendent random matrices, 

X = XiX2...X/, (50) 

where X^, i = 1,2,...,/, is rectangular of dimensions 
Ti X Ti_|_i, which tend to infinity in such a way that the 
rectangularity ratios 



■=■1 rr, 1 V"-^/ 

remain finite. I will follow the steps sketched in Sec. lIB 3l 
(cf. [ll|) to derive the mean densities of its singular values 
and if si = 1 also the eigenvalues. 



2. Eigenvalues of the product X assuming it is square and 
has rotationally-symmetric mean spectrum 

If si = 1, i.e., X is square, one may also ask about 
its eigenvalues. Assume that its mean spectrum has the 
rotational symmetry around zero (jl8p . Combining the 
A^-transform conjecture ([22)) with the multiplication law 
(1571) . one finds the appropriate formula. 



3. Eigenvalues of the product X assuming all X^ are square 
and have rotationally-symmetric mean spectra 



1. Singular values of the product X via cyclic shifts and 
multiplication law 

Let me for simplicity set / = 2 here. I begin from the 
eigenvalues of the Hcrmitian T3 x T3 random matrix. 



X^X = XTxjXiXa 



(52) 



However, as a first step, consider instead the T2 x T2 
matrix. 



Y = fxlXi) fXaX^) , 



(53) 



which differs from the previous one only by a cyclic shift 
of the terms. Therefore, as follows from ([T6l) . their A^- 
transforms are related by 



^Xtx(2) = ( 



(54) 



Now, Y is a product of two free Hermitian random matri- 
ces (their freeness is what I actually mean by the assumed 
"independence" of Xi and X2), and thus the multiplica- 
tion law ([25]) can be applied. 



N^{z) 



-Y^xtx,W^x.xtW = 
-Y^xtx.W^^xtx. 



(55) 



where in the last step I used an analogue of (|54p . Finally, 
inserting ([55]) into (l54l) , one obtains a multiplication law 
which allows to calculate the singular values of X from 
the singular values of Xi and X2, 

This formula may be generalized to any J, 



^xtx(^) 



[cf. Eq. (58) of the first position in [ll] 



(57) 



If moreover all = 1, i.e., all X^ are square, and also 
the mean spectra of all X^ are rotationally symmetric 
around zero, then the right-hand side of the multiplica- 
tion law (155)) may be expressed through the eigenvalues 
of Xi by virtue of the A^-transform conjecture ([22]) . which 
yields simply 



(59) 



9ix(z) = n^x,(^ 



The models T, W, V all possess the property ([T| 
since they are products of terms which separately exhibit 
this symmetry. Hence, one is allowed to use the multi- 
plication law ([5^ for them. Eventually, the rotational 
symmetry will be confirmed by Monte Carlo simulations. 



B. Product T 

1. Master equations for T 

Eigenvalues. It is now straightforward to write 
down the master equations for the nonholomorphic M- 
transform 9Jl = 9JlT(i?^) of the product T (gl)— they are 
comprised of the multiplication law ([5^ . 



^T(z)=n^s,(^), 

which links the J sets (|34ap - 

1=1 
_ z(z+l) 
'~ 9Ts,(z)' 
_ M,i {M,i + 1) 

— — ■ , I — 1, A, . . . , 1^ 

\w,l\ 



(60) 



(61a) 
(61b) 
(61c) 



for j — 1,2, . . . , J; these are (1 + 2 J + ^j) polyno- 

mial equations. They are valid inside the mean spectral 
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domain 2?, whose borderline (one or two centered circles) 
IS given by ([35a|) - (|35b)) . 

Singular values. These are obtained from the same 
master equations, albeit with O^t(-z) in replaced by 

^7VTtT(^) m- 

I will now simplify the above equations and solve them 
either analytically or numerically in a number of special 
cases, comparing the findings with Monte Carlo simula- 
tions. 



2. Example 



Eigenvalues. Before that, however, let us remark that 
if all the terms Sj have identical lengths, Lj = L, and se- 
quences of weights, Wji = wi, then (|61ap - (j61cp imply that 
all O^Sj (z) are equal to each other. This along with (pO)) 
in turn imply that the nonholomorphic M-transforms of 
T and any Sj are related by a simple rescaling of the 
argument. 



(62) 



Singular values. Unfortunately, there is no such scaling 
relation here; it is clear from trying to follow the above 
argument in conjunction with 



3. Example 1 

As the first example, take arbitrary J and Lj, but 
suppose that the weights for each Sj are equal to each 
other and denoted by 



for some J constants Wi 



I = 1,2, 



(63) 



Eigenvalues. Then, the master equations (|61ap - ()61cp 
give explicitly 

^s,(^) = |u;,f ^±3^, J -1,2,..., J, (64) 

which inserted into (j60p and after performing the func- 
tional inversion (f21]) leads to a polynomial equation of 
order J for 9Jt, 



i?2 



n^r(f + 



(65) 



In particular, if all Lj — L, Eq. ([55]) can be solved ex- 
plicitly, yielding (IT9)) . 



priR)-^2\w\'-j[l 



(I 



w 



2 _ fl^/J 



(67) 



for R < \w\, and zero otherwise. This could also be ob- 
tained by using the scaling relation (15^ and the proper 
result for S, i.e., Eq. (58) of ^. 

Singular values. The master equation for the holomor- 
phic Af -transform M = MxtT('Z) is a small modification 
of according to 



w 



(68) 



which is polynomial of order ( J + 1). 

The above findings, along with the erfc conjecture i 
are tested against Monte Carlo data in Fig. [T] [(a), (b), 
(c)]. 

Remark 1: Divergence at zero. Let me note that from 
the above equations one may easily derive an interest- 
ing feature of the level densities — their divergence close 
to zero. I will apply the logic presented in to this 
and the following models. Assume that for 2: — >■ 0, also 
z© ^ and zG ^ (recall, 971 = z© - 1, Af = zG - 1), 
which will be verified a posteriori. Then the denomi- 
nators of the right hand sides of [ (|65|) . (|68l) ] tend to 
nonzero constants, and consequently, © ~ z^/'^^^z^/'', 
G~z-^/(''+i), i.e.. 



pT{R)-R- 



i?-^ 0, 
X -> 0, 



where 



J. 



(69a) 
(69b) 

(70) 



[The initial suppositions thus hold true: z© ^ R^/'^ — > 
and zG ^ — > 0.] Note that the mean density of 
the singular values diverges for any J, while for the eigen- 
values, the radial mean spectral density is finite for J = 1 
(i.e., the model S, for which it zeroes) or J = 2, in which 
case however, the proper density ([5]) diverges. 

Remark 2: Is T really different from S ? Notice that 
if one opens up the brackets in the definition of the 
model T [([4|), ([1])] (let me focus on this Example 1 and 

all Wi — I), T = nf=i(-^7^^^ Sf^i Uj/), one obtains a 
• - - l/2\ 



where for short, w = Jl^^i w,. [Notice a certain simi- ^""^ °f 11^=1 Lj terms (multiplied by L i/2)^ each 



larity to the counterpart Eq. (IM)) for P; indeed, I will 
show (Sec. Ill C 21) that T of Example 1 is in some sense 
"inverse" to a certain class of models P.] 

Inserting ajt = or On = -1 (l35a |) - ((35b| ( no zero 
modes here, a = 0) into (j65p . one finds that the mean 
spectral domain is a centered disk of radius 



Ro^t. = \w\ 



(66) 



of which is a product of some J CUE matrices. Now, it 
is known that a product of CUE matrices still belongs 
to the CUE, hence, T looks like the model S of length 
L — except the fact that now the various terms are not 
statistically independent (free). So one might ask how 
relevant these correlations between the L terms are, i.e., 
how much the level densities of T and of S of length L 
differ. Figure [5] illustrates the theoretical mean spectral 
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FIG. 1: Theoretical level densities versus Monte Carlo data for the model T. Top row concerns the eigenvalues, middle row 
the eigenvalues plus the erfc form-factor, bottom row the singular values. Left column illustrates Example 1, right column 
Example 2. 



densities of these two random matrix models — they look 
completely different, showing that even though T may be 
recast as a sum of CUE matrices, correlations between 
them are important, and it is not our model S at all. 
Especially, the behavior at zero of the level densities in 
these two cases is entirely distinct, cf. Remark 1 above. 
Also, with growing L, the density of S tends to that of 
the square GinUE distribution (cf. [3), Ps'^ {R) 2R 
inside the unit circle and zero otherwise, while the limit 
of T is very different. 



4. Example 2 

As the second example, take arbitrary J, and let all 
the lengths Lj = 2, but with arbitrary weights Wji, Wj2- 
Eigenvalues. In this case, the master equations (|61ap - 
rcl) can be transformed into 



1 



1 



- 1 



{\w,i\-\wj2\r 



(71) 
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FIG. 2: Numerical verification of the property described in 
Remark 2 in Sec. Ill B 31 tliat even thougli tlie model T may 
be rewritten as a sum of CUE matrices, it is very different 
from tlie model S due to the correlations between the thus 
obtained CUE terms. 



which are quadratic for OTs (z)', they are supplied by the 
multiplication law (j60p . 

The borderline of the mean spectral domain is in gen- 
eral a centered annulus of radii (|35ap - (j35bl) . 



^ext. = n(k.i|' + k.2|'), (72a) 

(72b) 



J 



Rfnt. = n P^ll^ " I^i2| 

which reduces to a disk only if the (absolute values of 
the) two weights in at least one term are equal to each 
other. 

Singular values. As above, with the left-hand side of 
(|60)) changed according to (|22)) . 

These master equations are solved numerically and 
compared to the Monte Carlo eigenvalues in Fig. [T] [(d), 
(e), (f)]. 

Remark: Divergence at zero. Assume that there is a 
certain number 1 < W < J of terms Sj whose (abso- 
lute values of the) two weights are equal to each other, 
\wji \ = \wj2\; then one may ask about a behavior of the 
above level densities close to zero. For each of these 
W terms, Eq. ([7T|) becomes linear and yields (I already 
replace z by OK = z<5 — 1 ov M = zG — 1, respec- 
tively, and assume z@ — and zG — > for z — )■ 0), 
0^Sj(9H) = A\wj\'^z&/{zi& + 1), which tends to zero as zC5 
for z — > 0. For the remaining (J — W) terms in T, 
Eq. ([7Tt implies that any 9Tsj (2t) tends to a nonzero 
constant — |wj2p| for z — > 0. Substituting these 

results into , one discovers that the densities diverge 
at zero again according to (|69ap - (|69b|) . but with 

d = W. (73) 

(This finding is consistent with z0 — ^ and zG for 
z 0.) 



C. Product W 



Master equations for W 



Eigenvalues. The master equations for the mean spec- 
tral density of W (O are easily obtained through the mul- 
tiphcation law (15^1) . ^-wiz) — ^t{z)^p{z), from for- 
mula (p4)) (obviously with ri = 1) — they read 



K 



m-w{z) = cj^{z + 1) n ~ + M ^"r^-^) 



(74) 



where the rotationally-symmetric nonholomorphic A^- 
transform of T is given by (pn)) - (l61cp : one should sub- 
stitute here 971 = DJlw{R^) in the place of z. 

Singular values. Choose arbitrary ri , use the multipli- 
cation law ([55)) along with ([M)) — obtaining. 



wtw 



(75) 

where the holomorphic A^-transform of T^^T is given 
by [ (|22l) . (|60l) - (|61c|) ]: one should replace here z by 
M = Mwiwiz). 



2. Example 1 



Consider the situation described in Sec. Ill B 3l [i.e.. ar- 
bitrary J and Lj , weights for each Sj equal to each other 
(1631) ] plus arbitrary rectangularity ratios r^. 

Eigenvalues. The master equation is polynomial of or- 
der (J -I- a:), 



i?2 



(M+1) 



, nt2(f + i) 
n/^i(f + i) 

valid inside the disk of radius (|37ap - (|37b[) . 



i?nxt. = \w\a. 



(76) 



(77) 



[The internal radius vanishes because the density of the 
zero modes, a = 1 — minjrfc}, hence, either 9Jt = —1 or 
at least one term in the product over k in (ffS)) with 971 = 
a — 1 is zero. This same argument will be valid in all the 
models investigated henceforth.] 

Singular values. The master equation is polynomial of 
order (J-t- AT-f 1), 

J+i 



(M+i)(f + i) nt2(af + i 



M 



n^i( 



M 
riLj 



(78) 

Figures [3] [(a), (b), (c)] and|4][(a), (b), (c)] present nu- 
merical solutions to these equations and the correspond- 
ing Monte Carlo results, both in perfect agreement. 
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VIodel W; example 2 (Wjj+tVj2=l); singular values 




J=2(iv,,=0.2, H'2i=0.3) 

— J=3 (M.'n=0.2, H.2,=0.3, M.'3,=04) 

Z J=5 (wii=0.2, iV2i=0.3, W3i=0.4 
W4i=0.5, M'5i=0.6) 



FIG. 3: Theoretical level densities versus Monte Carlo data for the model W, with K — 2. Top row concerns the eigenvalues, 
middle row the eigenvalues plus the erfc form-factor, bottom row the singular values. Left column illustrates Example 1, right 
column Example 2. 



Remark 1: Divergence at zero. Following the same line 
of reasoning as in Remark 1 in Sec. lIIB~3l one recognizes 
that the divergences close to zero of the level densities 
stemming from [ ([75]) . ([75])] — disregarding their possible 
zero-mode part, which amounts to considering ri > 1, 
which means that there cannot happen Lj = 1/ri, i.e., 
one may disregard in our analysis the products over 
j in the above formulae — are governed by the number 
7^ = # {2 < A: < if : Tfe = 1}, and given again by (|69ap - 



(l69b| but with 

d^{J + l)5r,,i+'R. (79) 

Remark 2: T as an "inverse" of P with integer rect- 
angularity ratios. The above master equations of W im- 
ply the following peculiar joint property of the models 
T (of Example 1) and P: Assume that P (of length K; 
set for simplicity a = 1) is such that a certain number 
1 < J < K of the rectangularity ratios r^/ri are inte- 
gers greater than 1. If one multiplies this P from the left 
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FIG. 4: Theoretical level densities versus Monte Carlo data for the model W, with K — 5. Top row concerns the eigenvalues, 
middle row the eigenvalues plus the erfc form-factor, bottom row the singular values. Left column illustrates Example 1, right 
column Example 2. Figure (f) does not include theoretical graphs due to numerical complications encountered in solving the 
master equations. 



by a product T of length J such that the lengths Lj of 
its terms are equal to these integer rect angularity ratios 
(and Wj = 1), then the master equations [ ([75]) . (|78p ] turn 
into 



i?2 = (an + 1) 



.7+1 



n 

2<k<K: 
Tk is not integer > 1 



— + 1 

rk 



(80) 



M 



n 



2<k<K: 
is not integer > 1 



.7+1 



M 



(81) 



Comparing [ (|M)) . ([5T|) ] with (|24p . one recognizes that 
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R X 



FIG. 5: Numerical verification of tfie property described in Remark 2 in Sec. Ill C 21 Tfie tfieoretical grapfis are tlie level 
densities of the model Pi, obtained from (I24|) with all rt = 1- 




X 



FIG. 6: Theoretical level densities [the eigenvalues (a), the eigenvalues plus the erfc form-factor (b), the singular values (c)] 
versus Monte Carlo data for the model V. 



they are the master equations for a model P = P1P2, matrix T [of length {K — 1) and the lengths of its terms 

where Pi = Bi . . . Bj, with Bj being the square GinlJE Lk — rk/ri] yields 
random matrices ©, while P2 is the model P with the 
terms Afe for which r^/ri — Lj removed from it. 

In particular, if the model P is such that all of its rect- 
angularity ratios r^/ri, fc > 2, are integers greater than j.ad 2 2/K-1 

1, then multiplying it from the left by the appropriate 9Jl = i? — 1, i.e., (i?) = — i? , (82) 
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(Af + l)(M + l) 

M 



K 



(83) 



In other words, the model T (of Example 1) is "m- 
verse" (at least on the level of the mean densities of 
eigenvalues and singular values) to the model P with 
rectangularity ratios r^/ri, k > 2, equal to the lengths 
Lj, with the "unity" being Pi (US])]. These results 

are verified numerically in Figs. [S] [(a), (b)]. 



3. Example 2 



In the setting of Sec. Ill B 41 [arbitrary J, all Lj — 2, 
arbitrary weights Wji, Wj2\ with arbitrary r^, the master 
equations l(7i)) and ([75l) should be supplied by the result 
(|7ip . These can be solved numerically and are compared 
to the numerical simulations in Figs. [3] [(d), (e), (f)] and|3] 
[(d), (e), (f)]. 

An analytical result can be obtained for the mean spec- 
tral domain — it is a disk of radius (|37ap - p7bl) . 



^n(i 



Wji\ 



\Wj2 



(84) 



This should be contrasted with the analogous situation 
for the model T, where the domain is in general an an- 
nulus ([72a|) - (f72bl) . 

Remark: Divergence at zero. Combining the reason- 
ings from Sees. Ill B 41 and III C 21 one finds that for the 
mean spectral domains to touch zero, there must be a 
number W of weights obeying \ wji \ = \wj2\, in which case 
the divergences at zero of the level densities (without tak- 
ing into account possible zero modes) are still described 
by ((69a|) - ((69bl) but with 



d^{yV +l)5r,,l+'R. 



(85) 



Notice that this is consistent with ([801) . and even sug- 
gests how an expression for d might look like for a most 
arbitrary model W. 



Wi = TiP, of the form of Example 1 above (Sec. Ill C 2|) 
[i.e., the lengths Jt of arbitrary, with arbitrary lengths 
Lij of Sij but with the weights Wiji = Wijj ^fLij in- 
dependent of l\ the lengths Ki of P^, variances af^ 
and rectangularity ratios r^fc of A^fc also arbitrary], 
having moreover arbitrary rectangularity ratios (jSlI) . 
Si = N^i/Ni^K,+i = nf'=i''»'i (of course, Si = 1 if one 
investigates the eigenvalues of V) . 

Eigenvalues. The master equation is polynomial of or- 
der + 



|w|2a2 



n 



m 

Si 



Ji+1 



i+lTik 



1 



(86) 



where a = OLi nf=2 ^*fe and w = OLi 11^1 and 
it is valid inside the disk of radius p7a|) - (|37bp . 



(87) 



Singular values. The master equation is polynomial of 
order ELliJ^ + K^) + I, 

M + 1 



Figures[n][(a), (b), (c)] illustrate numerical solutions to 
these polynomial equations and test them against Monte 
Carlo simulations, finding perfect agreement. 

Remark 1: Generalized Bures products. Notice that 
for / = 1, Eqs. [(ESI), (Ell)] reduce to Also, 
setting / = 1, Ji = 1, Ki = l,i/7 = l,(T = l,si=:ri = l, 
ill = 2 in (EH) yields a cubic equation, whose solution is 
the original Bures distribution (|1^ . Of course, an even 
broader generalization of the Bures model would be to 
include distinct weights Wiji in the S^j; I refrain from 
explicitly writing down the pertinent master equations, 
even though it is a straightforward step from the above 
results. 

Remark 2: Divergence at zero. The level densities 
stemming from [ ((5S)) . (EH)] diverge at zero (beyond possi- 
ble zero modes) again according to (|69ap - (|69bp but with 



D. Product V 



Master equations for V 



The mean densities of the eigenvalues or singular values 
of any product V ^ can be directly derived from (f74|) or 
(|75l) using the multiplication laws (|58|) [or (|59|) ] or ([57]) . 



c^-I]((^. + l)'5s.,l + 7^.), 



(89) 



where for short, TZi = {2 < k < Ki : rik — 1/si+i}. 
Formula (EH) reduces to ((75)) for / = 1, as it should. 



III. CONCLUSIONS 



2. Example 



A. Summary 



I will consider in this paper just one instance of 
the model V, namely, any length / but all the terms 



In this paper, I attempted to show how free probability 
theory (the multiplication law in the realm of Hermitian 
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random matrices and the addition law for non-Hermitian 
ones) greatly simplifies otherwise difficult calculations of 
the mean densities of eigenvalues and singular values of 
the generalized Bures products T, W, V, which are ran- 
dom matrix models of relevance in quantum information 
theory. 

B. Open problems 

As already mentioned, this article is certainly only the 
initial step (nevertheless important, especially from the 
point of view of quantum information theory) in learning 
about the models T, W, V. A major endeavor would be 
to consider finite matrix dimensions and attempt a com- 
putation of the complete JPDF of the eigenvalues of these 
models. One could also investigate some of their univer- 
sal properties (cf. Sec. II A3")) . Actually, inspired by p^ . 
one could check whether the level densities themselves 



show any sign of universality, just like for P. 

A pressing challenge is to prove the three hypotheses 
which my derivation is founded upon: the iV-transform 
conjecture (cf. Sec. II B 2p . the single ring conjecture 
(cf. Sec. II B 5)) . and the erfc conjecture (cf. Sec. II B 6)) . 

It would also be desirable to understand more about 
the features of our models important for quantum entan- 
glement theory, e.g. analyze their von Neumann entropy 

(ED. 
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